Dynamics of Quintessence Models of Dark Energy with Exponential Coupling to the 

Dark Matter 



o 
o 

(N 
00 



> 

(N 
(N 

o 
o 



Of 
i 1 

o 



X 



Tame Gonzalez, Genly Leon and Israel Quiros 1 

1 Universidad Central de Las Villas, Santa Clara CP 54830, Cuba 
(Dated: February 5, 2008) 

We explore quintessence models of dark energy which exhibit non-minimal coupling between the 
dark matter and the dark energy components of the cosmic fluid. The kind of coupling chosen 
is inspired in scalar-tensor theories of gravity. We impose a suitable dynamics of the expansion 
allowing to derive exact Friedmann- Robertson- Walker solutions once the coupling function is given 
as input. Self-interaction potentials of single and double exponential types emerge as result of our 
choice of the coupling function. The stability and existence of the solutions is discussed in some 
detail. Although, in general, models with appropriated interaction between the components of the 
cosmic mixture are useful to handle the coincidence problem, in the present study the coincidence 
can not be evaded due to the choice of the solution generating ansatz. 
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I. INTRODUCTION 

At the time, all the available observational evidence 
suggest that the energy density of the universe today is 
dominated by a component with negative pressure of, al- 
most, the same absolute value as its energy density. This 
mysterious component, that violates the strong energy 
condition (SEC), drives a present accelerated stage of 
the cosmic evolution and is called "dark energy" (DE). 
A pleiad of models have been investigated to account for 
this SEC violating source of gravity. Among them, one 
of the most successful models is a slowly varying scalar 
field, called 'quintessence' 0,0. 

Many models of quintessence assume that the back- 
ground and the dark energy evolve independently, so 
their natural generalization are models with non-minimal 
coupling between both components. Although experi- 
mental tests in the solar system impose severe restric- 
tions on the possibility of non-minimal coupling between 
the DE and ordinary matter fluids @, due to the un- 
known nature of the dark matter (DM), it is possible 
to have additional (non gravitational) interactions be- 
tween the DE and the DM components, without conflict 
with the experimental data. 1 Since, models with non- 
minimal coupling imply interaction (exchange of energy) 
among the DM and the DE, these models provide new 
qualitative features for the coincidence problem 0, @] . It 
has been shown, in particular, that a suitable coupling, 
can produce scaling solutions that are free of the coinci- 
dence. The way in which the coupling is approached is 
not unique. In reference [91, for instance, the coupling 
is introduced by hand. In [fj the type of coupling is not 
specified from the beginning. Instead, the form of the in- 
teraction term is fixed by the requirement that the ratio 



1 It should be pointed out, however, that when the stability of 
dark energy potentials in quintessence models is considered, the 
coupling dark matter-dark energy is troublesome 



of the energy densities of DM and quintessence had a sta- 
ble (non-zero) equilibrium point that solves the problem 
of the coincidence. In 0], a suitable interaction between 
the quintessence field and DM leads to a transition from 
the matter domination era to an accelerated expansion 
epoch in the model of Ref.[(|. A model derived from 
the Dilaton is studied in [1]. In this model the coupling- 
function is chosen as a Fourier expansion around some 
minimum of the scalar field. 

A variety of self-interaction potentials have been stud- 
ied in DE models to account for the desired evolution that 
fits the observational evidence. Among them, a single 
exponential is the simplest case. This type of potential 
leads to two possible late-time attractors in the presence 
of a barotropic fluid @, [Io|: i) a scaling regime where the 
scalar field mimics the dynamics of the background fluid, 
i. e., the ratio between both DM and quintessence energy 
densities is a constant and ii) an attractor solution dom- 
inated by the scalar field. Given that single exponential 
potentials can lead to one of the above scaling solutions, 
then it should follow that a combination of exponentials 
should allow for a scenario where the universe can evolve 
through a radiation-matter regime (attractor i)) and, at 
some recent epoch, evolve into the scalar field dominated 
regime (attractor ii)). Models with single and double ex- 
ponential potentials has been studied also in references 

Hp. 

Aim of the present paper is to investigate models with 
non-minimal coupling among the components of the cos- 
mic fluid: the dust dark matter and the quintessence 
field (a scalar field model of dark energy). To specify 
the kind of coupling, we take as a Lagrangian model, 
a scalar-tensor (ST) theory with Lagrangian written in 
Einstein frame variables. Otherwise, one might also con- 
sider the Lagrangian with additional non-gravitational 
coupling between the matter species as an effective the- 
ory. 

In this paper, to derive solutions of the field equations 
of the model, we use an ansatz that makes possible easy 
handling of the differential equations involved at the cost, 
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however, of loosing the possibility to avoid the coinci- 
dence problem through an appropriate choice of the inter- 
action between the DM and the DE (quintessence field). 
Actually, we impose the dynamics of the expansion by ex- 
ploring a linear relationship between the Hubble param- 
eter and the time derivative of the scalar field. 2 Using 
this relationship we can solve the field equations explic- 
itly. However, as it will be shown, unlike other models 
where the coincidence problem is solved (or smoothed 
out) through the choice of an appropriated interaction 
between the DM and the quintessence field, in the present 
investigation, solutions where the ratio between the dark 
matter and the quintessence energy densities is a constant 
~ 1 (it is the necessary requirement to avoid the coinci- 
dence problem), are unstable so that the coincidence do 
arises. 

Here we consider a flat Friedmann-Robertson- Walker 
(FRW) universe filled with a mixture of two interacting 
fluids: a background of DM and the quintessence field. 
Since there are suggestive arguments showing that ob- 
servational bounds on the "fifth" force do not necessar- 
ily imply decoupling of baryons from quintessence [l3j . 
then baryons are to be considered also as part of the 
background DM interacting with the quintessence field. 
In other words, like in reference [l3j we are considering 
universal coupling of the quintessence field to all sorts 
of matter (radiation is excluded). Since the arguments 
given in the appendix of reference [TJj to explain this pos- 
sibility are also applicable in the cases of interest in the 
present study, we refer the interested reader to that refer- 
ence to look for the details. However we want to mention 
the basic arguments given therein: A possible explana- 
tion is through the "longitudinal coupling" approach to 
inhomogeneous perturbations of the model. The longi- 
tudinal coupling involves energy transfer between matter 
and quintessence with no momentum transfer to matter, 
so that no anomalous acceleration arises. In consequence, 
this choice is not affected by observational bounds on 
"fifth" force exerted on the baryons. Other generaliza- 
tions of the given approach could be considered that do 
involve an anomalous acceleration of the background due 
to its coupling to quintessence. However, due to the uni- 
versal nature of the coupling, it could not be detected by 
differential acceleration experiments. Another argument 
given in [l3| is that, since the coupling chosen is of phe- 
nomenological nature and its validity is restricted to cos- 
mological scales (it depends on magnitudes that are only 
well defined in that setting) , the form of the coupling at 
smaller scales remains unspecified. The requirements for 
the different couplings that could have a manifestation 
at these scales are that a) they give the same averaged 
coupling at cosmological scales, and b) they meet the 



2 This relationship is the simplest possible and, in the absence 
of any other information is, in a sense, the most natural since 
the Hubble parameter fixes the time scale. This argument was 
suggested to us by Diego Pavon. 



observational bounds from the local experiments. We 
complete these arguments by noting that they are ap- 
plicable even if the coupling is not of phenomenological 
origin, like in the present investigation where the kind of 
coupling chosen is originated in a scalar-tensor theory of 
gravity. 

The paper has been organized as follow. In section 
II the details of the model are given. The method used 
to derive FRW solutions is explained in section III. We 
use the 'two fluids' approach: a background fluid of DM 
and a self-interacting scalar field (quintessence). Spe- 
cific exponential couplings with different exponents, that 
are inspired in ST theory and lead, correspondingly, to 
self-interaction potentials of single exponential and dou- 
ble exponential class, are studied separately. In section 
IV a study of existence and stability of the solutions is 
presented. We point out that, due to the choice of the 
solution generating ansatz, the coincidence problem can 
not be evaded since, a constant ratio between the en- 
ergy densities of the components of the cosmic mixture 
is never a stable attractor. In section V, conclusions are 
given. 



II. THE MODEL 

We consider the following action that is inspired in a 
scalar-tensor theory written in the Einstein frame, where 
the quintessence (scalar) field is coupled to the matter 
degrees of freedom: 

s = [ d^x^MilR-hvctf -V(4>) 

+C 2 {4>)C (matter) }. (1) 

In this equation R is the curvature scalar, <p is the 
scalar (quintessence) field, V{4>) - the quintessence self- 
interaction potential, C 2 (0) - the quintessence- matter 
coupling function, and £ ma tter is the Lagrangian density 
for the matter degrees of freedom. This action could be 
considered, instead, as an effective theory, implying ad- 
ditional non-gravitational interaction between the com- 
ponents of the cosmic fluid. 

Although in the present study we will be concerned 
mainly with FRW spacetimes with flat spatial sections, 
here for generality we write the FRW line element in the 
form: 



ds 2 = -di 2 + a(t) 2 ( — T + rW), 

1 — kr l 

di! 2 = dO 2 + sin 2 9dcf> 2 , (2) 

where t is the cosmic time, r, 9, <p are the spatial (radial 
and angular) coordinates and k is the spatial curvature, 
which we take to be zero in this investigation. We use 
the system of units in which 8nG = c = h = 1. 
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The spacetime is filled with a background pressure- 
less dark matter fluid and a quintessence field (the scalar 
field <j)). As already said the baryons (a subdominant 
component at present, but important in the past of the 
cosmic evolution) are included in the background of dark 
matter. In the introductory part of this paper we have 
already commented on the possibility of a universal cou- 
pling of the dark energy to all sorts of matter, including 
the baryons (and excluding radiation). 

The field equations that are derived from the action 
© are: 

31Z 2 + ™ = Pm + U 2 + V, (3) 
a z 2 

2H + 'SH 2 + \ = (1 - 7 )p m - \ft + V, (4) 
a 2 

+ 3H0 = -V' + (lnX)'p m , (5) 

where we have introduced the reduced notation X{(f)) = 
C(^)(37-4)/2 The parameter 

7 is the barotropic index 
of the background fluid (DM) . The " continuity" equation 
for the background is: 

p m + 3jHp m = -(In X)'<j> p m , (6) 

or, after integration 

p m = Ma-^X-\ (7) 

where M is a constant of integration. In the former equa- 
tions the dot accounts for derivative in respect to the 
co-moving time t, while the prime denotes derivative in 
respect to <f>. We now proceed to derive exact solutions 
to the above field equations by fixing the dynamics of the 
expansion. 

III. DERIVING SOLUTIONS 



corresponding differential equation. As it will be shown, 
exponential coupling functions with different exponents 
represent the simplest situations to study in the present 
model. Correspondingly, the ansatz © implies that only 
self-interaction potentials of the exponential form (in- 
cluding their combination) can be considered. Integra- 
tion of equation ([5]) implies that 

a = (9) 

where the scale factor has been normalized so that the 
integration constant has been absorbed into it. After 
this, the equation ([7]) can be rewritten in the following 
form: 

p m {(t>) = Me-^X- 1 ^). (10) 
If one adds up equations ([3]) and ^ then one obtains: 

H + 3H 2 = ^ Pm + V. (11) 

Substituting ([5]) in (0) , and comparing the resulting equa- 
tion with (fTTj) . a differential equation relating V and the 
coupling function X (and their derivatives in respect to 
the scalar field variable (j>) can be obtained: 

dV 1Tr ,„ , JN fdlaX A(7-2)\ 

Consider further equation (|8|) written in the form d(f> — 
Ae?(ln a) . It is then worthwhile rewriting of the equation 
(fl"2f in the following form: 

v + *v-(£-^)£*-. (13) 

where now the prime denotes derivative in respect to the 
variable N — In a and we have substituted p m (<fi) from 
equation p0[) . Equation §5§ can then be integrated in 
quadratures: 



In order to derive exact analytic (flat; k = 0) solutions, 
one should either fix the dynamics of the cosmic evolution 
or fix the functional form of the self-interaction potential 
V{4>). In the present study we fix the dynamics of the 
model by imposing the following constraint 



(8) 



that involves the Hubble parameter and the square root 
of the kinetic energy of the quintessence field (A is an 
arbitrary constant parameter) . As it will be immediately 
shown, this relationship enables one to reduce the sys- 
tem of differential equations ([3][5]) to a single first order 
differential equation, involving the self-interaction poten- 
tial V together with its derivative in respect to the scalar 
field variable <p an d the corresponding derivative of the 
coupling function. Therefore if one chooses further the 
form of the coupling function, the functional form of the 
self-interaction potential can be found by solving of the 



dxp 



^Me-T^A-i^) + V{<f>) 



2A 2 



6- A : 



■(t + to),(14) 



or, if one introduces the time variable dr 



e 2 



i^X^^dt, 



M ■ 



L *X(cl>)V(<P) 



2A 2 



6 - A 2 



(r + ro). (15) 



In consequence, once the function X(<fi) (or X(a)) is 
given as input, then one can solve equation (|12p (or (| 13() ) 
to find the functional form of the potential V(4>) (or 
V(a)). The integral (HH) (or dT5j)) can then be taken 
explicitly to obtain t = t{4>) (or r = t(0)). By inversion 
we can obtain cf> = <fi(t) (or (f> = 4>{t)) so, the scale factor 
can be given as function of either the cosmic time t or 
the time variable r through Eqn.®. 
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A. Particular Cases 

We shall study separately the simplest situations that 
can be considered once the choice ([5]) is made. In both 
cases one deals with coupling functions of the exponential 
form: X(4>) = Xq exp (n<f>), where X$ and n are constant 
parameters. If one chooses n = (2 - 3^/2)/ y/u + 3/2 
where u) is the Brans-Dicke coupling parameter, then the 
action fl} corresponds to Brans-Dicke theory written in 
the Einstein frame (EF). In this case the EF scalar field 
<\> is related to the Jordan frame scalar field <f> through: 
d<j) = d(f>/ (ui + 3/2). In general, once the dynamics |0) is 
imposed, this kind of coupling function leads to double 
and single exponential potentials depending on r. The 
importance of this class of potentials in cosmology has 
been already outlined in the introductory part of this 
paper. 



1. Case A 

Let us consider the simplest situation when, in Eq. 
H; 



d\nX _ A(7- 2) 
dcj) ~ 2 : 



(16) 



i. e., we are faced with an exponential coupling function 
of the form mentioned above with r = \(-f — 2)/2 in the 
exponent. In this case the equation (|12[) simplifies to 



d ^ + xv = o, 



(17) 



which can be easily integrated to yield to a single expo- 
nential potential: 

V = V o e~ A0 . (18) 
In consequence the equation (fT5j) can be written as 



dw 



\/w 2 



A 2 



= (J,(t + T ) 



(19) 



where w — exp(- 



A 2 = V X /M and [i 



7^/M(A 2 — 6)/8, so we have the explicit solution 

4>{t) = 0o + In (sinh [/x(r + r )]^^) ) , (20) 
and, consequently: 

a{r) = a sinh [/x(r + r )] . (21) 

The dimensionless density parameter (for the i-th com- 
ponent fli = pi/'SH 2 ) and the Hubble expansion param- 
eter can be given as functions of the redshift z also. Ac- 
tually, if one considers that a(z) — ao/(l + z), where 



ao = a(z = 0) (for simplicity of the calculations we 
choose the normalization ao — 1) then; 



6 - A 2 (1 + z )37+7(6-A 2 )/2-A 2 

(1 + z )3 7 +7(6-A2)/2-A2 + A 2 ' 



(22) 



and 



H(z) = bJ^{1 + z)37+7(6-a^)/2 + (1 + (23) 



where B = y/2V Q /{& - A 2 ). Note that, Q m is a max- 
imum at z = oo: f2 m (oo) = (6 — X 2 )/6A 2 . In gen- 
eral one can write that, at the epoch of nucleosynthe- 
sis Q m (oo) = (1 — e), where e is a very small number 
(the small fraction of dark energy component during nu- 
cleosynthesis epoch) e = [6(A 2 - 1) + \ 2 }/6A 2 . Tak- 
ing into account the observational fact that, according to 
model-independent analysis of SNIa dataflj], at present, 
(z = 0); fi m (0) = 1/3, then (f22"|) can be rewritten as 



(1 - e)(l + z )37+7(6-A 2 )/2-A 2 
" m(z) = (1 + z )37+7(6-A2)/2-A2 + ( 2 _3e)' 

where now the solution exhibits only two free parameters 
A and e (the DM EOS parameter 7 = 1). 

Other physical magnitudes of observational interest are 
the quintessence equation of state (EOS) parameter and 
the deceleration parameter, that are given by the follow- 
ing expressions: 



W0 



-1 + 



q = 



= -1 



2 2 



(25) 



(26) 



respectively. 



2. Case B 

A second very simple choice is to consider, in Eq. (|12p : 
d 



(In X (</>)) = const = a X((f>) = X e°" t ', (27) 



where Xq is an integration constant. In this case, by 
integration of (fT2")) one obtains 



V((t>) = V Q e~ x * + WV~ (a+37/A)0 , 
where the constant 

M f\(2- 1 )-2a 



Wn = 



2X \a + 37/A - A 



(28) 



(29) 
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FIG. 1: The evolution of Q m (thick solid line) and fl^ (dashed 
line) vs z is shown for the model with a single exponential 
potential. The following values of the free parameters e = 
0.01 and A = 0.3 have been chosen. Equality of matter and 
quintessence energy density occurs approximately at z as 0.3— 
0.4. 



FIG. 3: We plot the evolution of the deceleration parameter 
9 vs z for the Model A (single exponential potential). As 
before, the following values of the free parameters e = 0.01 
and A = 0.3 have been chosen. Again, the curve with the 
smallest value of the parameter A is the one that fits better 
in the observational data set on a present stage of accelerated 
expansion. 
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FIG. 2: We plot the dynamical EOS parameter of the 
quintessence vs z, for the Model A (single exponential po- 
tential), for three different values of the parameter A: 0.3 
(thick solid line), 1.41 (solid line), and 2.24 (dashed line) re- 
spectively. In all cases the free parameter e is chosen to be 
e = 0.01. Note that only the curve with the smallest value 
of the parameter A (thick solid line) meets the requirements 
of the present observational data favoring a value of the EOS 
parameter u)^ ~ — 1. 



We are faced with a self-interaction potential that is a 
combination of exponentials with different constants in 
the exponent. The usefulness of this kind of potential 
has been already explained in the introduction and will 
be briefly discussed later within the frame of the present 
model, when we study the stability of the corresponding 
solution. If one introduces the variable y — e~2<^ the 
equation (|15p can be written as 



dy 



Vv 2 + 1> 2 



2 A 2 XqVq 
6 - A 2 b 2 



(r + r ) (30) 



where b 2 = (X Vq)/(M + W X Q ) and I = a + 3 7 /A - A. 
Integration of ((30)) yields to: 



4>{t) = <f) + In (sinh [/z(t + to)]) 
and, consequently: 

a(r) = a (sinh [fj,(r + tq)])~ 2 ' 1 . 



(31) 



(32) 



where fi = — | J g 2 ^ X °Y° ■ On the other hand, by using 

the Friedmann equation ([3]) and inserting the potential 
(|28p. one obtains the Hubble parameter as function of 
the quintessence field <j>: 



2V 
6 - A 2 



b 2 e -(a+3~ ( /\)4> + e -A0 



(33) 



where we have considered Eqn. ([T0|) . In terms of the scale 
factor one has, instead, 



H\a) = 



2V 
6- A 2 



(34) 



Using the redshift variable z the above magnitudes can 
be written as follows: 



H 2 (z) 



2V b 2 
^A 2 



(^ + 1) 



(aA+3 7 ) 



+ 



(l+^) > 
b 2 



(35) 



In turn, the energy density of DM can be written as 
Pm (z) = (M/X )(z + 1) qA+3 t = Po (z + 1)^+37, so 
the DM dimensionless density parameter can be given 
as function of z also: 



(^) 



(6 - A 2 )p (z + l)«A+3 7 -A 2 

2V b 2 (z + l)«A+3 7 -A2 + l/ b 2 ■ 



(36) 
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In order to constrain the parameter space of the so- 
lution one should consider the model to fit the obser- 
vational evidence on a universe with an early matter 
dominated period and a former transition to dark en- 
ergy dominance. 3 Therefore A 2 — aX — 37 < 0, besides, 
as before, one can write f2 m (oo) = (1 — e), where e is a 
small number (the small fraction of dark energy compo- 
nent during nucleosynthesis epoch). In correspondence 
b 2 = [(6 — A 2 )/(l — e)](p /2V ). Besides, if one considers, 
as before, that according to model-independent analysis 
of SNIa dataQ, at present (z = 0); fl m (0) = 1/3, then 
1/b 2 = 2 — 3e. After this Eqn. (f3"o]) can be rewritten in 
the following way: 



n m (z) = (1 - e) 



(z + l)«^+37-A 2 



(z + l)« A + 3 ^ 2 + (2-3e)' 



(37) 



so the solution depends on three free parameters: a, A 
and e (the barotropic index of the DM is fixed: 7 = 1). 
These can be chosen so that the solution fits well the ob- 
servational data. Other physical magnitudes of observa- 
tional interest are the quintessence EOS parameter and 
the deceleration parameter in equations ([25)) and ([26]) . 
respectively. 

In Figure [TJ we show the evolution of both dimen- 
sionless DM and scalar-field energy densities f2 m and 
respectively vs z for the single exponential poten- 
tial (Case A). The corresponding evolution for Case B is 
very similar. In Figure^ the evolution of the dynamical 
quintessence EOS parameter vs z is plotted for the Model 
A (single exponential potential), for three different val- 
ues of the parameter A: 0.3 (thick solid line), 1.41 (solid 
line), and 2.24 (dashed line) respectively. In all cases the 
free parameter e is chosen to be e = 0.01. Note that 
only the curve with the smallest value of the parameter 
A (thick solid line) meets the requirements of the present 
observational data favoring a value of the EOS parameter 
ujfj, ~ —1. Meanwhile, in Figure O we plot the evolution 
of the deceleration parameter q vs z for the Model A. 
Again, the curve with the smallest value of the param- 
eter A is the one that fits better the observational data 
set on a present stage of accelerated expansion [l7l li~8j . 

Figures Q] and [5] show the behavior of = u^>(z) and 
q = q(z) for different values of the free parameter a: 
0.1 (thick solid line), 1 (solid line), and 5 (dashed line). 
The free parameters e and A have been fixed ad hoc, 
guided by the results of the study of case A (e = 0.01, 
A = 0.3). It is apparent that the present values of 
and of q, do not depend on a, so that this parameter can 
not be determined from the observational data on the 
present stage of accelerated expansion of the universe. 



3 Since we are dealing with models with non-minimal coupling 
between the quintessence and the matter fields, here we refer, 
mainly, to a model-independent analysis of SNIa observational 
datarig.ri5l.ril. 
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FIG. 4: We plot the dynamical EOS parameter of the 
quintessence vs z, for the Model B (double exponential po- 
tential), for three different values of the parameter a: 0.1 
(thick solid line), 1 (solid line), and 5 (dashed line) respec- 
tively. In all cases the remaining free parameters are chosen 
to be e = 0.01 and A = 0.3 respectively. At small redshift the 
dependence upon a is only weak. 
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FIG. 5: We plot the evolution of the deceleration parameter 
q vs z for the Model B, for three different values of the pa- 
rameter q: 0.1 (thick solid line), 1 (solid line), and 5 (dashed 
line) respectively. As before, the following values of the free 
parameters e = 0.01 and A = 0.3 have been chosen. The 
present value of the deceleration parameter is independent on 
the value of a. 



This parameter could be of impact if the model is applied 
to study early stages in the cosmic evolution. 



IV. EXISTENCE AND STABILITY OF THE 
SOLUTIONS 

We now turn to the study of the stability of the so- 
lutions found. In order to keep the study as general as 
possible, we do not specify any concrete model for the 
dark energy. We only require that the DE EOS parame- 
ter u>de > — 1- For this end we rewrite the field equations 
in the following form. The Friedmann equation (we 
include the most general situation with spatial curvature 
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X 


V 


Existence 


Stability 


1 


1 


always 


stable node if < A < saddle if v2 < A < \/6; unstable node otherwise 


X 2 

fi 


1 


< A < VE 


unstable node 


A^ 

2 


A^ 

2 


always 


saddle point if < A < y/2; stable node if A > \J~2 



TABLE I: The properties of critical points for case A. 



Point (x, y) 


Ai 


A 2 


(1,1) 


A^-6 


A 2 - 2 


K -§- L =x 


3 — i— 
6 , 2 


A^ + 2 

2 


(hi hi) 

\ 2 ' 2 / 


A 2 +2 
2 


2 — A 2 



TABLE II: The eigenvalues for Case A. 



X 


2/ 


Existence 


Stability 


1 


1 


always 


stable node if < A < \/2 and < a < 3 ~ A A ; saddle point if either < A < \/2 and 
a < 3 ~^ A , or \/2 < A < \/3 and < a < 3 ~^ A ; unstable node otherwise 


|(A + a) 


1 


< A < V3 
and < a < 2^ 


saddle point if < A < \/2 and i < a < 
unstable node otherwise 


2 


A 2 
2 


always 


stable node if A > \/2 and < a < j-\ unstable node if < A < \pl and 
a > i-; saddle point otherwise 



TABLE III: The properties of critical points for case B. 



Point (x, y) 


Ai 


A 2 


(1,1) 


A 2 - 2 


A 2 + Aa - 3 


(A(A + a)/3,l) 


3 — A 2 — Aa 


1 - Aa 


(hi hi) 

\ 2 ' 2 > 


Aa- 1 


2- A 2 



TABLE IV: The eigenvalues for Case B. 



3H 2 + 3^r = p r , 



(38) 



the Raychaudhuri equation: 

2H - 2-^ = -{p m + p m + V4> + P4> ), (39) 

the continuity equation for the quintessence field: 

P4 , + 3iJ( j o + pf) = -Q, (40) 

and the continuity equation for the background matter: 

p m + 3H(p m +p m ) = Q, (41) 

where the dot accounts for derivative with respect to the 
cosmic time and Q = — (In X)'p m is the interaction term. 
Through this section the prime will denote derivative in 
respect to the new variable TV = In a, that is related with 
the cosmic time through dN = Hdt. 



Let us introduce the following dimensionless phase 
space variables: 4 



x = fi^, ; y = Sltot = fi<£ + 



(42) 



After this, the Friedmann equation (|38[) can be rewritten 
in the following way: 



k 



a 2 H 2 ' 



(43) 



The governing equations ([315 |) .([3"9" 1) .(|4"0" ]) . and (gTJ) can 
be written in terms of the phase space variables in the 
following way: 

x = (y - x)(lnx)' + x(y - 1 + 3cj (.t - 1)) 

V' = 0/- l)(y + 3^x). (44) 



4 See reference [Tgj ] for an alternative treatment of a related sta- 
bility study. 
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The above equations represent an autonomous system of 
equations if (ln%)' and uj^ do not depend on N explicitly. 
In the remaining part of this section, for simplicity, we as- 
sume that it is the case, so that the system ([4"4"]l is an au- 
tonomous one. Besides we consider x( a ) — Xo & S , where 
S is some constant parameter. This choice of the cou- 
pling function comprises many useful situations (includ- 
ing the solutions we have derived before) and it implies 
that (lnx)' = 5. We assume, also, that uj m = ^> 7 = 1, 
i. e., the background fluid is dust. For the flat space case 
(k = 0), the system (|44p should be complemented with 
the following constraint equation: 



0.8- 



0.6- 



0.4 ■ 



0.2- 



0.2 



0< y-x < 1, 



(45) 



which follows from requiring that the positive dimension- 
less matter energy density parameter VL rn < 1. 

The first step towards the study of the dynamics of 
the autonomous system (|44|) is to find its critical points 
(x c ,y c ) (x',y') — (0,0). Then one can investigate 
their stability by expanding equations (|44|) in the vicinity 
of the critical points x — x c + u, y = y c + v (up to terms 
linear in the perturbations u,v): 



(46) 



where A is the matrix of the coefficients in the expan- 
sion. The general solution for the evolution of the linear 
perturbations can be written as: 




u = Ul e x > N +u 2 e X - N , 
v = vi e AlJV + v 2 e x * N , 



(47) 



where Ai and A 2 are the eigenvalues of the matrix A. 

In table [J we show the properties of the critical points 
(including existence and stability) for the case A above, 
meanwhile, in table [III the corresponding eigenvalues are 
given. In this case X — Xq e~ x ^/ 2 . As seen, for the 
values of A that make the model fit better the obser- 
vational data set (0 < A < 1), the first critical point 
(x,y) — (1,1) => f2^ = 1 (quintessence dominated 
phase) is a stable attractor, while the scaling solution 
(x,y) = (A 2 /6, 1) => fWft = 6/A 2 - 1 > 1, is always 
unstable. The phase portrait for this case is shown in 
figlE All of the trajectories in phase plane (x, y), diverge 
from the unstable point (matter dominated scaling solu- 
tion) and converge towards the attractor (quintessence 
dominated) solution. 

The properties of the critical points and the corre- 
sponding eigenvalues for case B, where X = Xo e a ^, are 
shown in tables |TTT] and |IV] respectively. It is apparent 
that, for the relevant ranges of the free parameters < 
A < 1 and 0<a<3/A— A, the quintessence dominated 
solution (first critical point (x, y) — (1, 1)) is always a sta- 
ble node. The (matter dominated) scaling solution- sec- 



FIG. 6: The phase plane for Case A (7 = 1, A = 0.3). The 
critical point (1, 1) is stable (a sink) so that the quintessence 
dominated solution is the late time attractor. The scaling 
regime (0.015, 1) is an unstable point. The saddle is located 
at (0.045, 0.045). All the phase space trajectories diverge 
from the unstable point and converge towards the attractor. 



ond critical point ((x,y) = (A(A + a)/3, 1))- could be ei- 
ther asaddle (1/A < a < 3/A— A) or , otherwise, an unsta- 
ble node. The phase portrait figJ7]shows that, for the val- 
ues of the free parameters chosen (A = 0.3, a — 5.7), all 
phase space trajectories diverge from the unstable node 
(third critical point (x,y) = (A 2 /2, A 2 /2)-quintessence 
dominated solution with curvature) and, either converge 
towards the stable attractor solution dominated by the 
quintessence (first critical point), or are repelled by the 
saddle point (second critical point- the matter domi- 
nated scaling solution). This result, that is generic for 
both solutions A and B, shows that, since the scaling 
(quintessence dominated) solution with f2 m /fi^ < 1 is 
not even a critical point of the corresponding autonomous 
dynamical system (equations l|44p). then the coincidence 
problem could arise in the cases studied in the present 
investigation. In the next subsection we will show in a 
more definitive manner that this is indeed the case. 



A. The coincidence problem 

Let us investigate whether, in the situations of interest 
in the present study, the question: why are the energy 
densities of the dark matter and of the dark energy of the 
same order precisely at present? arises. For this pur pos e 
it is recommended to study the dynamics of the ratio [lj]: 



P<i> 



7^ 



(48) 



in respect to the variable N = In a, that, as said before, is 
related with the cosmic time t through dN = Hdt. The 
following generic evolution equation holds for r: 



f(r), 



(49) 
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X 



FIG. 7: The phase plane for Case B (7 = 1, A = 0.3 and 
a = 5.7). The quintessence dominated solution (point (1, 
1) in the phase plane) is an attractor. The scaling solution 
(critical point (0.6, 1)) is a saddle, meanwhile (0.045, 0.045) 
is an unstable point. For different values of the parameter 
a, the separation between the saddle (scaling phase) and the 
stable (quintessence dominated) critical points varies. 



where the prime denotes derivative in respect to the vari- 
able ./V and / is an arbitrary function (at least of class 
C 1 ) of r. One is then primarily interested in the equi- 
librium points of equation (|49p . i.e., those points r e i at 
which f(r e i) — 0. After that one expands / in the 
neighborhood of each equilibrium point; r = r e i + e,, 
so that, up to terms linear in the perturbations e,: 
f(r) = {df/dr^^+Oiti) => e'i = {df/dr)^. This 
last equation can be integrated to yield the evolution of 
the perturbations: 

ei = e 0i exp[(df/dr) rei N], (50) 

where eoi are arbitrary integration constants. It is seen 
from (|50p that, only those perturbations for which: 

(df/dr) rei < 0, (51) 

decay with the time variable N, and the corresponding 
equilibrium point is stable. The coincidence is evaded if 
the point p m / P<t> — f e i < 1 is stable. 

Now, if we take into account the conservation equa- 
tions (|4Tj|) . and (|4Tj) , where Q = — (lnX)'p m , and since 
Q<P = + 1), implying that (see equation (|2l)|) ) the 
quintessence EOS parameter uj^ can be written as func- 
tion of r: ujcf, — — l + A 2 (r + l)/3, then, for the cases under 
study here, the function / can be given by the following 
expression: 

/(r)=r{(A 2 -<5)(r + l)-3}, (52) 

where 6 = (lnX)' = nX. For the case A: S = — A 2 /2, 
while for the case B: 5 — Xa. It can be easily checked 
that, in both cases the only stable equilibrium point is the 



one for which the ratio r e o = 0, i.e., it is the quintessence 
dominated solution, so that the coincidence can not be 
evaded. 



V. CONCLUSIONS 

We have found a new parametric class of exact cos- 
mological scaling solutions in a theory with general non- 
minimal coupling between the components of the cos- 
mic mixture: cold dark matter and dark energy (the 
quintessence field). Baryons, although subdominant at 
present, in the past of the cosmic evolution played a 
relevant role and could be included as part of the dark 
matter component in the present setup. There are sug- 
gestive arguments showing that observational bounds on 
the "fifth" force do not necessarily imply decoupling of 
baryons from quintessence [l3T] . To specify the general 
form of the coupling we considered a scalar tensor theory 
of gravity written in the Einstein frame. An alternative 
interpretation is to consider it as an effective theory, im- 
plying additional non-gravitational interaction between 
the components of the cosmic fluid. 

In order to derive exact flat FRW solutions to the 
model under study, we have assumed a linear relation- 
ship between the Hubble expansion parameter and the 
time derivative of the scalar field. Mathematically this 
assumption allows to reduce the original system of sec- 
ond order differential equations to a single, first order 
differential equation, involving only the self-interaction 
potential and the coupling function (together with their 
first derivatives in respect to the scalar field variable). 
However, simplicity of the mathematical handling is at 
the cost of retaining the problem of the coincidence. In 
fact, the solution generating ansatz (jHJl, implies that solu- 
tions where the ratio p m j p^ = const ~ 1 are not stable. 
In the cases studied, only the quintessence dominated 
solution is stable. Anyway, the assumed relationship be- 
tween the square root of the scalar field kinetic energy 
and the Hubble parameter is the simplest possible and, 
in the absence of any other information is, in a sense, the 
most natural since the Hubble parameter fixes the time 
scale. 

We have concentrated our study in exponential class 
of coupling functions. Brans-Dicke theory is a particular 
member in this class. Exponential coupling functions can 
lead to self-interaction potentials of the following class: 
A) single exponential potential and B) double exponen- 
tial potential. The stability and existence of the solutions 
found have been also studied. For this purpose we have 
applied a fairly general method in which one does not 
need to specify any model for the dark energy. In both 
cases (case A and case B), the dynamical system exhibits 
three critical points. For the values of the free parameters 
that are allowed by the observational data, the scalar field 
(quintessence) dominated solution is always an attractor, 
meanwhile the scaling (matter dominated) solution can 
be either an unstable node or a saddle point. 
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We conclude that models with non-minimal coupling 
between the dark energy and the dark matter are easy 
to handle mathematically if one assumes a suitable dy- 
namics. This is some times, at the cost of retaining the 
problem of the coincidence. The present investigation 
could be complemented by the study of different classes 
of self-interaction potentials and/or the choice of suitable 
dynamics allowing to evade the coincidence problem; one 
of the motivations to study models with interaction be- 



tween the components of the cosmic mixture. 
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